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ABSTRACT 



Binary sequences find increasing use in electrical 
engineering applications of ranging, time measurement and 
communications. A property of interest in these applications 
is the autocorrelation function of the binary sequence or 
pair of sequences. Of the 2^ possible sequences of length n, 
only a few have usable autocorrelation functions. There is, 
to date, no procedure known which will provide the sequence 
having a specific autocorrelation function, except for very 
particular cases . 

In this report, known properties of complementary sequences 
are reviewed. Almost complementary sequences are defined and 
the procedure to obtain them is outlined. A formula is 
derived for the number of different autocorrelation functions 
of the 2^ possible sequences of length n bits. A computer 
search is implem.ented with the objective of discovering se- 
quences with desirable autocorrelation functions. 
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I. INTRODUCTION 



This study is concerned with binary sequences and their 
autocorrelation functions. The objective is to obtain auto- 
correlation functions with sidelobe levels less than or equal 
to predetermined values . Sequences or groups of sequences 
which can provide this property are very attractive for use 
in systems whose performance depends on the autocorrelation 
function magnitude. Such systems are used in communications, 
ranging and time measurement. The problem is to find these 
sequences or "good" codes and the rules to construct them, if 
such rules exist,. 

A class of such codes are the complementary sequences. 

They are pairs of sequences with the characteristic that the 
sum of their autocorrelation functions is a waveform that has 
no sidelobes . These sequences were first considered by 
M. J. Golay [Ref. l] and further investigated by S. Jauregui 
[Ref. 2]. They are used here to develop a technique for con- 
structing another class of "good" codes, the almost com.ple- 
mentary sequences. 

A. PLAN OF THE RESEARCH 

The efforts to solve the problem follow two directions: 

(1) Experimentation with the binary sequences and 
their properties, to find the rules which give desirable auto- 
correlation functions . 
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(2) Computer search of binary sequences of several 
lengths to find the ones with small autocorrelation sidelobes. 

As a result the following were achieved: 

(1) Discovery of the almost complementary sequences. 
These are pairs of binary sequences which can be constructed 
using complementary sequences . Their autocorrelation functions 
when added have sidelobes of predetermined magnitude, polarity 
and position. 

(2) Discovery of codes of lengths 5^20 which have auto- 
correlation functions with sidelobes less or equal to one. 

For example, for length n = 20, only three such codes were 
found and for n = 15, none. The computer programs used here 
can be used to select codes with any sidelobe levels. 

Other results of interest are: 

(1) Derivation of the formula for the number of the 
different autocorrelation functions of all sequences of 
length n. 

(2) Construction of computer programs which can be 
used in other cases as well. For example, an algorithm for 
the automatic production of all possible binary numbers of 
length n was devised. This algorithm can be used to select 
codes having certain properties, such as codes with a fixed 
number of ones and zeros. 

B. PLAN OF THE REPORT 

Chapter II provides the necessary background by giving 
the definitions and basic properties of the binary sequence 
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and autocorrelation function. The formation of the autocorre- 



lation function R^(t) of a two level function v(t) is developed. 
Matched filters and their realizations are also discussed. 

Complementary sequences and their basic properties are 
reviewed in Chapter III, to form the basis for the material on 
almost complementary sequences . 

In Chapter IV, the almost complementary sequences are 
defined, and rules for their construction are given. 

The formula for the number of different autocorrelation 
functions in all sequences of length n is developed in 
Chapter V . 

Computer programs and their algorithms are discussed next 
in Chapter VI . 

Chapter VII gives a few applications. 



13 



II. BACKGROUND 



In this section some basic concepts are discussed and 
definitions given as background material. 

A. BINARY SEQUENCE 

A binary sequence is a list of elements each of which can 
have one of two distinct values. These values are usually 
represented either by +1 and -1 or by 1 and 0. Often when 
the +1 and -1 convention is used, the ones are omitted and 
only + and - are written. 

For example, sequence A can be written: 

A = ++-+ or A = 1101 or A = +1+1-1+1 

The number of elements in a sequence is the length denoted 
here by n. For the above example, n = 4. 

The voltage equivalent v(t) of a binary sequence is a 
time waveform where 1 is represented by a voltage level +V 
and 0 by a voltage level -V . 

For the sequence A = 1101, v(t) is given in Fig. 1 where 
£ is the bit duration. 
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v(t) 



+V 



1 10 1 



£ 2e 3£ 4e 



t 



-V 



Fig. 1. Voltage equivalent of a binary sequence. 

In this study the 1,0 notation is used. Also, when 
sequence is used, binary sequence is implied. 

3. AUTOCORRELATION FUNCTION 

The autocorrelation function of a two-level, time-limited 
voltage v(x) is defined as the integral 



R^^(t) is a measure of the similarity between a voltage or 
signal and its phase shifted version where all values of time 
delay t are considered. 

The way to find the autocorrelation function R^(x) of a 
digital sequence v(t), is to "slide" the sequence past itself 
to the right and left and at each position form the product 
of the sequence and its shifted replica. Then the area of 




/ 



oo 



v( t ) v( t-T ) dt . 
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the product waveform is taken and this gives the autocorrelation 
of the sequence at this position. 

The procedure is illustrated in Fig. 2, by showing the 
complete steps for two "shift" positions t = 0 and x = s. 

It can be seen from the equation for R^(t) and from Fig. 2 
that when v(t) is a piece-wise constant function, R^(t) will 
be a piece-wise linear. The linear segments terminate at 
multiples of e, a bit duration. 

The autocorrelation is an even function which has its 

maximum value at x = 0. gives the level of the main 

lobe. Secondary maxima are the sidelobe levels. For the 

2 

example of Fig. 2, the mainlobe level is e. 

Sidelobe levels are R^(£) = -V^e and R^(3e) = +V^e . 

Actually the shape of the autocorrelation function of a 
sequence v(t) is obtained easier by letting v = 1 and e = 1. 

For example, to form the autocorrelation function of the 
sequence A = 1101, the sequence is written and its delayed 
version is placed beneath. For example, x = 0 gives 

I 

1101 ; 

1101 ’ 

In each position the elements of these two similar se- 
quences are compared. If they are the same (both zeros or 
ones) they form an agreement; if not (one zero, one one) they 
form a disagreement. The number of disagreements is subtracted 
from the number of agreements and the result is proportional 
to the autocorrelation function at this position. Here there 
are only four agreements and the result is 4. 
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V 



2 



v( t) v( t ) 



e 2e 3e 4e t 

(a) T = 0, R^(0) = 4V^£ 






< 

rt 


/( t-e ] 


1 


-v^ 


e T 


£ 4 


£ t 



(b) X = e, R^(e) = -V^e 




Fig. 2. Autocorrelation of a binary sequence. 
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Now a shift is made as follows 



T 

f 



1101 ! 
110 '1 



and by the same method there are one agreement and two dis- 
agreements, so the result is -1. 

Similarly, the next position is 



1101 

11 



01 



which gives 0. 

The next position 



1101 

1 



101 



gives 1. 

And finally. 



the 



last po 



T 

llOlj 

t 

t 



sition 



1101 



always gives 0 . 

Since the autocorrelation is an even function, R^(t) = 

If T is considered to provide a shift to the future, 
-T is a shift to the past . 

It is not necessary for the values given by shifting to 
the left to be written, because they are the same with the 
ones resulting from shifting to the right. So, by convention 
the autocorrelation of the sequence A can be written 



R^(t) = 4, -1,0, 1,0. 

This convention will be followed in the rest of this report. 
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C. MATCHED FILTERS 



Some interesting properties of matched filters will be 
listed here. These properties are derived in the literature 
[Ref. 3,4]. 

A matched filter is the best linear filter for detection 
of a pulse signal v(t) in noise. The impulse response h(t) 
of such a filter is a delayed, time inverted replica of the 
input. If v(t) is the input to the matched filter, its impulse 
response is 

h(t) = Mv[-(t-to)], t^O 
h(t) = 0 , t< 0 . 



where M = an arbitrary constant 



t = time delay inherent in the filter, 
o 



From linear system theory, the output g(t) of a matched 
filter is 

r +00 



:(t) = M 



v(t-T)v [-(T-t^)]dX 
o 



Let X = t-T to obtain 



g(t) = m/ v(X ) v(-t+X+to)dX 

\J —00 

r 

or g(t) = M I 



v(X)v(X- 5 )dX = M 



where C = t-to. 

The output is maximized when 5 = 0 ot t = t^. 

It is concluded then, that the output g(t) of the matched 
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filter is the autocorrelation function of the input. Fig. 3 
illustrates the concept. 



g(t) = M R^(t-to) 

Fig. 3. Matched filter. 

Matched filters for two-level voltages (binary sequences) 
can be realized by tapped delay lines or shift registers as 
shown in Fig. 4. The tapped delay line realization uses in- 
verters at these positions where a zero element occurs in the 
sequence. The shift register realization uses a reference 
register where the original sequence is stored. This can be 
a read only memory (ROM) for example. Another register 
receives the input sequence v(t). 

In both realizations, +1 units of current flow through 
the load resistor R, for each element of the input sequence 
that agrees with the "stored" sequence. And -1 unit of current 
flows through R, for each element of the input sequence that 
disagrees with the "stored" sequence. The net output current 
through R^ (and voltage across Rj^) is proportional to the 
number of elements which agree less the number of elements 
which disagree. The output g(t) is, therefore, a measure of 
the autocorrelation function of the input sequence. 

Since the systems of Fig. 4 perform a discrete comparison 
and summing instead of multiplication and integration, then 
g(t) is a discrete version of R^(t). For example, the sequence 
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g( t ) a n 




(b) Shift register 



Fig. 4. Matched filter realizations. 
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1101 has the autocorrelation function shown in Fig. 2, whereas 
the output g^(t) of the discrete matched filter corresponding 
to 1101 is as shown in Fig. 5. 



u 


V 


(t) 










1 




_4e -3s -2s 






2s 3 s 4 S 


“1 







Fig. 5. Discrete matched filter output. 

For a sequence of n elements, the peak output is n units 
of voltage. It is clear then, how signal detectability 
improves as n increases . 

All binary sequences have autocorrelation functions with 
sidelobes of various values . The sequences of interest are 
the ones with either small or negative sidelobes. The 
problem is to find these sequences. 

The next section addresses the issue of forming a pair of 
sequences or codes, which when processed with matched filters 
and the outputs added, yield a waveform with one mainlobe and 
no sidelobes . This scheme provides good detectability of 
binary sequences . 
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Ill . COMPLEMENTARY SEQUENCES 



This section reviews complementary sequences and their 
basic properties . 

A set of complementary series is defined as a pair of 
equally long, finite binary sequences which have the property 
that the number of pairs of like elements with any given 
separation in one series is equal to the number of pairs of 
unlike elements with the same separation in the other series . 
For example the two series: 

A = 00010010 
B = 00011101 

are complementary. In A there are three like elements (denoted 
by t below) separated by one element. 

Q^O^O 1 0^0 1 0 

In 3 there are; 

00011101 
a a a 

three unlike elements (denoted by a below) separated by one 
element . 

Similarly for all possible separations the number of like 
elements in A and unlike elements in B are as follows : 
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Separation 



Number of 
Likes in A 



Number of 
Unlikes in B 



1 

2 

3 

4 

5 

6 
7 



1 



1 



4 



2 



2 



3 

3 



3 

3 

4 
2 
2 



1 



1 



Series having the complementary property were conceived by 
Marcel J. E. Golay in connection with the optical problem of 
infrared multislit spectrometry [Ref. l^S]. 

Complementary series have interesting autocorrelation 
functions . If the autocorrelation of each sequence is taken 
and these two autocorrelations summed, the result is zero for 
all X except x = 0 . At x = 0 , the sum is twice that of either 
sequence. Therefore, the sum of the autocorrelation functions 
has one main lobe and no sidelobes. 

For example, consider the sequences or codes 



A = 00010010 



B = 00011101 



Sequence A has autocorrelation function: 



R^(x) = 8, -1, 0, 3, 0, 1, 0, 1, 0. 



Sequence B has autocorrelation function: 



Rg(x) = 8, 1, 0, -3, 0, -1, 0, -1, 0. 
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Their sum is: 



Z = + RgCt) = 16, 0, 0, 0, 0, 0, 0, 0, 0. 

This property can be treated in equation form as follows: 

Let a^ and (i = 1, 2, , n) be the elements of 

two sequences A and B each of length n. Assume a^ and b^ can 
be either +1 or -1. Then the respective values of the auto- 
correlation functions will be 



n-j 

c. = Z a. a. . for j^O. 
J i = l 



d . = E b . b . , . for j ^ 0 . 

^ i=l ^ 



Also cj 


= c j for 


A 

O 


d. 


= d . for 


j< 0 . 


3 


-3 



The necessary and sufficient condition for the pair of 
sequences to be complementary is : 



c . + d . 
] 3 



and 



+ d. 
3 



0 for j 0 
2n for j = 0 . 



where j ranges from -n+1 to n-1. 
Or in expanded form. 



n-j 

E 

i = i 



a • 



1 



i+j 



n-j 



+ E 
i-1 




0 for j 0 
2n for j = 0 . 
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When the elements of the sequence are 0 or 1 , then the 
autocorrelation function is obtained by modulo-two addition. 

In this case the necessary and sufficient condition for the 
series to be complementary is that 

n-j n-j 

I (a^© = Z (b^ @ b^^j © 1) for all j, 1< j ^ n-1 

i=l i=l 

This complementary property can be tested experimentally 
by using matched filters since the output of a filter matched 
to its input is the autocorrelation function of that input. 

This realization is shown in Figure 6 . The output of the 
system of Figure 6 is shown in Figure 7 for the complementary 
sequences 

A = 111-1 

B = 11-11 

A. GENERAL PROPERTIES 

1 . Number of Elements 

The number of elements in two complementary series 
are equal. If they were different, the pair of extreme elements 
of the longest series would remain unmatched by an unlike pair 
of elements with the same spacing in the other series. 

2 . Symme try 

Two complementary series (A, B) are interchangeable (B,A); 
that is, one can take the place of the other. This results from 
the symmetry of the definition with respect to two complementary 
series . 
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ig. 6. Scheme for testing the complementary property. 
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Fig. 7. Analog matched filter processing of 
two complementary sequences. 
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Sequence Length Even 



3 . 



A necessary condition for a sequence pair to be com- 
plementary is that their length n be an even number. 

4 . Sequence Length Sum of Two Squares 

Another necessary condition for a pair of sequences 
to be complementary is that their length be the sum of the 
squares of two integers . The proof was developed by 
S. Jauregui [Ref. 4]: 

5 . Transformations 

A single pair of complementary series can be rhe basis 
for the construction of 64 pairs of complementary series. 

a. Order of Complementary Sequences 

Denote the reverse of A by A^ . For example, 
if A = 1110, then A^ = 0111. The order of the elements of 
either or both of a pair of complementary series may be 
reversed. This follows from the fact that by reversing a 
sequence its autocorrelation function remains the same . The 
proof is developed in Section V. 

b. Complementing the Sequence 

Denote the complement of A by A. For example, if 
A = 1110, then A = 0001. One or both of a pair of complementary 
sequences can be complemented — putting zeros in the place of 
ones and ones in the place of zeros, without affecting their 
complementary property. This follows from the fact that by 
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complementing a sequence, its autocorrelation function remains 
the same. Section V provides the proof. 

c. Complementing Elements of Even Order 

Denote the complement of the even order elements 
of A by Ag . For example, if A = 1110, then A^ = 1011. Com- 
plementing the elements of even order in each sequence — putting 
zeros in the place of ones and ones in the place of zeros, does 
not affect their complementary property. 

It is concluded from the above properties that a single 

pair (A,B) of complementary sequences can be the basis for the 

0 

construction of 2° = 64 pairs of complementary series (some of 
which might be identical) by either performing or not performing 
the following six operations: 



a . 


Interchanging 


the sequences . 


b. 


Reversing the 


first sequence. 


c . 


Reversing the 


second sequence. 


d. 


Complementing 


the first sequence. 


e . 


Complementing 


the second sequence . 


f . 


Complementing 
sequence . 


the elements of even order of each 



As an example, consider the complementary pair A = 00010010 
and B = 00011101: 

Applying a gives B = 00011101 and A = 00010010. 

Aptlying b gives A^ = 01001000 and B = 00011101. 
Applying c gives A = 00010010 and B^ = 10111000. 
Applying d gives A = 11101101 and B = 00011101. 
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Applying e gives A = 00010010 and B = 11100010. 

Applying f gives A^ = 01000111 and B = 01001000. 

6 e 

By applying the above properties properly, the original 
pair can be reproduced: 

B = 00011101 and A = 00010010 

B^ = 10111000 and A = 00010010 

B^ = 10111000 and A^ = 01001000 

3^ = 01000111 and A^ = 01001000 

(f^)e = A = 00010010 and (A^)^ = B = 00011101. 

The last pair is the same as the original one . 

6 . Allowable Lengths 

Since, as was mentioned before, the number of elements 
in complementary sequences must be even and equal to the sum 
of two squares, the allowable sequence lengths up to 50 are 

2, 4, 3, 10, 16, 18, 20, 26, 32, 34, 36, 40, SO. 

It has been verified by trial, though, that complementary 
sequences for length 18 do not exist. 

7 . Hamming Distance 

The Hamming distance of two binary sequences A and B 
is defined as the number of positions in which these two 
binary sequences differ. This can be written in modulo two 
notation as follows: 
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For 



i=n 

D(A,B) = Z a. © b. 
i = l ^ 

example, the Hamming distance of the two sequences 

A = 0100 
B = 1111 

i = 4 

D(A,B) =Z a. ©b. =(1+0+1+ 1) =3 
i = l ^ 

Now for a complementary pair of sequences, it has been 
proven that their Hamming distance is always = y [Ref. 2]. 

For example, the complementary pair of length n = 10 

A = 1001010001 
B = 1000000110 

has Hamming distance 

D(A,B) =1=5. 

8 . Kernels 

A Kernel is a basic length sequence which cannot be 
decomposed into shorter length sequences . The shortest 
possible complementary pair is 11 and 10. This pair or any 
of each transformation, which was mentioned before, is called 
a kernel of length two or a quad. 
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Some possible Kernel lengths are 
2, 10, 18, 26, 34, 50. 

It might be the case, though, that complementary pairs for 
some of them do not exist. For example, lengths n = 4 and 
n = 8 have complementary pairs, but are not Kernels because 
they can be constructed from n = 2 and n = 4 sequences 
respectively . 

Among all the above mentioned Kernel lengths, it has 
been verified by M. J. Golay that n = 18 does not exist. 

Also, it has been verified by S. Jauregui [Ref. 2] 
through exhaustive computer search that for n = 26 only the 
Kernel shown in Table I exists, not taking into account all 
allowable transformations. 

For n = 34 a non-exhaustive computer search by 
S. Jauregui revealed no Kernel. An exhaustive computer search 
was not possible, due to the great computer tim.e required. 



It is possible, however, that a complete search for 



n = 34 


could be achieved 


in the future. 


using 


new techniques . 


The following table shows the Kernels of 


n = 


2, 10, 26, 


ignoring allowable trans 


formations . 










Table I 






n 


Kernels 
Number of 
Kernels 


of Length 2, 10, 
A Sequence 


26 


B Sequence 


2 


1 


10 




11 


10 


2 


1001010001 

0101000011 




1000000110 

0000100110 


26 


1 


010011011110101111 

00111010 


101100100001111111 

00111010 
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Note: (a) The possible Kernel of n = 18 does not exist. 

(b) Partial computer search for n = 34 found no 
Kernel . 

9 . Number of Ones in Complementary Sequences 

S. Jauregui showed [Ref. 2] that the equation 
n = (n-p-q)^ + (p-q)^ 

holds for two complementary sequences A and B of length n 
where p is the number of ones in A and q the number of ones 
in B . 

This leads to the conclusion that the number of ones 
in each of the sequences of a complementary pair cannot be 
arbitrary, but has to satisfy the above relation. 

For example, for complementary sequences of length 
n = 2 the number of ones in A and B can be respectively 

either (2, 1) 
or ( 1 , 0 ) . 

B. SYNTHESIS 

If the sequences A,B are complementary, they can be used 
to generate other complementary pairs as follows: 

(a) If A = a^ a 2 a 3 a^_^ a^ 

B = b^ b2 bg ^n-1 ^n 

are a complementary sequence pair, then the sequences 



34 



*' ■ ®1 ^2 % ’^2 

D = aj — -a^ bj -— 

are also complementary. 

For example, consider the complementary pair: 

A = 0001 with autocorrelation = ‘+5 0, - 1 , 0 

B = 0010 with autocorrelation RgC^) = 4, -1, 0, 1, 0. 

Then the sequences 

C = 00010010 
D = 00011101 

are also complementary with autocorrelation functions 
Rj^(t) = 8 , -1, 0, 3, 0, 1, 0, 1, 0 

R^(t) = 8 , 1, 0, -3, 0, -1, 0, -1, 0. 



(b) If A = 


3 ., 3 « 


a , a 




12 3 


n -1 n 


B = 




b , b 




X 2 3 


n -1 n 


are a complementary 


sequence 


pair, then the sequences 



C = a. 


b-. 


I-) b 


b 


1 


1 


2 2 


n n 


D = a. 


b. 


a 0 bo 


3 


1 


1 


2 2 


n n 



are also complementary. 

For example, consider the same sequences 
A = 0001 
B = 0010 
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The sequences 



C = 00000110 
D = 01010011 

are also complementary with autocorrelation functions 

R(.(x) = 8, 3, 0, 1, 0, -1, 0, 1, 0 
Rj^(t) = 8, -3, 0, -1, 0, 1, 0, -1, 0. 

(c) If (A,B) (C,D) are two complementary sequences 
pairs, A of length n and C of length n, then the pair 

= a'^I A^2 A^n B^l B^2 B^n 

V2 = A^n A^l B^n B^l 

is also complementary, where if an exponent is one the A or B 
sequence is left unchanged, whereas if the exponent is zero 
the A or B sequence is complemented. 

For example, consider the complementary pairs 

A = 11 
B = 10 

and 

C = 00 
D = 01. 

Then the pair 

= 00000110 
V 2 = 11001010 
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M 




is also complementary with autocorrelation functions 

R (t) = 8, 3, 0, 1, 0, -1, 0, 1, 0. 

^1 

R (t) = 8, -3, 0, -1, 0, 1, 0, -1, 0. 

^^2 



The above methods make possible the generation of comple- 
mentary • sequence pairs of greater lengths than the original 
ones. They can be applied in succession to generate very long 
sequences which are very useful in many applications. 

For example such a complementary pair used in a communica- 
tions system with matched filter processing like that of Fig. 6 
can improve signal detectability in the presence of considerable 
noise since the summer output voltage will consist of a large 
main lobe and no sidelobes. 



C. SUPPLEMENTARY AND CYCLIC COMPLEMENTARY SEQUENCES 

Complementary sequences are subsets of two larger sets, 
namely supplementary and cyclic complementary sequences . 

1 . Supplementary Sequences 

Consider two sequences 



A = ai a2 a3 a^_^ a^ 



B 



b 



n-1 



b 

n 
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Now let 



^ -- 



— a 



n-1 



II = a a ^ 
n n-2 



a 



III = b3 b3 



2 

b 



n-1 



IV = b b . 
n n-2 



b. 



The expression of the sequence pair A, B in the form 
(I, II, III, IV) is called sequence quadruple. Supplementary 
sequences are quadruples of sequences with the property that 
the total number of likes at each spacing equals the total 
number of unlikes at the same spacing. 

In terms of their autocorrelation function, the sum 
of the four autocorrelation functions is zero any place but 
T = 0 , where it is four times the length of the sequences . 

For example. 



A = 1001010001 
B = 1000000110 

Writing in (I, II, III, IV) form gives 

I = 10000 

II = 10110 

III = 10001 

IV = 01000. 
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The new sequences (I, II, III, IV) have auirocorr’e- 
lation functions 

Rj(x) = 5, 2, 1, 0, -1, 0 
Rjj(x) - 5, —2, -1, 2, -1, 0 
Rjjj(t) - 5, 0, -1, -2, 1, 0 
Rj^^(t) = 5, 0, 1, 0, 1, 0. 

The sum of the autocorrelation functions is 

E = Rj(x) + + ^iii^’’-^ 

= 20 , 0 , 0 , 0 , 0 , 0 . 

In this example the sequences A, B are complementary. 
In general they do not have to be though, since the supplementary 
sequences are a larger set. This is illustrated in the following 
example . 

Consider 

I = 000100111 
II = 000101001 
III = 000101000 
IV = 000110110 

with autocorrelation functions 

Rj(x) = 9, 2, -1, 0, 1, 0, -3, -2, -1, 0 

Rjj(x) = 9, -2, 1, 0, -1, 2, 1, 0, -1, 0 

= 9, 0, 3, -2, 1, 0, 3, 2, 1, 0 

Rj^(x) = 9, 0, -3, 2, -1, -2, -1, 0, 1, 0. 
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The sum is 



E = Rj.(t) t + RjjjCt) + Rjy(t) 

= 36, 0, 0, 0, 0, 0, 0, 0, 0, 0 

which is zero except at the position x = 0, where it is 36. 

From the sequences (I, II, III, IV) the sequences 
A, B can be constructed 

A = 010000110001101010 
B = 000101100111000000 



In this example neither the A, 3 sequences nor 
the (I, II, III, IV) ones are complementary, which demonstrates 
the fact that supplementary sequences are a larger set. 

Figure 8 gives another example of the supplementary 

property . 



2 . Cyclic Complementary Sequences 



In general a cyclic sequence is a never ending periodic 
sequence of zeros and ones which has period of n elements. A 
cyclic complementary sequence pair is a pair of cyclic 
sequences, each of period n, where the number of likes in one 
sequence equals the number of unlikes in the other one, for 
all possible spacings . In equation form 

n n 



C. 

3 




© a 



n- j +1 



= E b ® b . © 1, l^j^n-1. 

i + 1 J ^ 



In terms of autocorrelation functions the sum of the 
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"two periodic autocorrelation functions is zero except at 
X = kn, where k = 0, 1, 2, .... 

For example, the cyclic complementary sequences 

A = 01010011, 01010011, 

B = 00000110, 00000110, 

have periodic autocorrelation functions 

R^(t) = 8, -4, 0, 0, 0, 0, 0, -4, 8, .... 

Rg(x) = 8, 4, 0, 0, 0, 0, 0, 4, 8, .... 

The sum is 

Z = R^(t) + RgCt) = 16, 0, 0, 0, 0, 0, 0, 0, 16... 

Figure 9 demonstrates the above example. 

Complementary sequences are always supplementary and 
cyclic complementary, but the opposite is not always true. 
Supplementary and cyclic complementary sequences constitute 
a larger set. 
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IV. ALMOST COMPLEMENTARY SEQUENCES 



Complementary sequences are attractive for use in communi- 
cations and ranging systems because the usable receiver output 
voltage has no sidelobes. Sequences providing "small” side- 
lobes compared to the main lobe may also be useful in some 
applications . 

It may even be desirable to have a small sidelobe at a 
known position and level. For example, such a sidelobe can 
be used to measure doppler as shown in Section VII. 

Two binary sequences whose summed autocorrelation functions 
exhibit two small sidelobes are called almost complementary in 
this report, since they exhibit properties similar to comple- 
mentary sequences. 

It should be noted that there are many sequences with small 
sidelobes, but the ones of interest here, are those for which 
certain rules hold. By applying these rules, the magnitude 
and position as well as polarity (positive or negative) of the 
sidelobe can be precisely predicted. By knowing these rules 
almost complementary sequences or codes can be constructed. 

An extensive search for "good" almost complementary 
sequences was made. As a result, it is possible to list general 
rules for constructing almost complementary sequences with pre- 
dictable sidelobe levels, polarity and position. There is a 
distinction between positive sidelobe and negative sidelobe 
sequences . 



1+4 



A. POSITIVE SIDELOBE SEQUENCES 



Here two cases are developed. 
+ 1 

1 . Sidelobes ~ K Away 



If A, B are complementary sequences of length n and 
the new sequences C and D of length K = 2n are formed as 
follows , 



C = AA 
D = BB 

the sum of the autocorrelation functions of C and D gives only 
two positive sidelobes of magnitude half that of the main lobe 
at T = - j K. 

For example, if n = 8, then K = 16 with the main lobe 
level of 16, and sidelobe levels of 8 at t = - 4. 

For example, consider the two complementary sequences 
of length n = 4, 



A = 0001 
B = 1011, 

and construct the almost complementary sequences 

C = AA = 00010001 
D = BB = 10111011 

having autocorrelation functions 
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R^(t) = 8, 1, 0, -1, 4, 1, 0, -1, 0 
Rj^(t) = 8, -1, 0, 1, 4, -1, 0, 1, 0. 

Adding gives 

S = Rq(x) + Rq(t) = 16, 0, 0, 0, 8, 0, 0, 0, 0. 

So there are only two positive sidelobes with levels half the 

+ 1 + 

level of the main lobe and at a distance x = -- 2 K = -4. 

Figure 10 illustrates the above example. 

Another example is 

A = 0001001011100010 
B = 0100011110110111 

A, B are complementary of length n = 16. Constructing C = AA, 

D = BB gives 

C = 00010010111000100001001011100010 
B = 01000111101101110100011110110111 

C, B are almost complementary of length K = 32 with auto- 
correlations 

Rq (t) - 32, — 1, 0, — 1, 0, 3, 0, — 5, 0, — 1, 0, 9, 0, 
1, 0, 1, 16, -1, 0, -1, 0, -1, 0, -1, 0, -3, 
0, 5, 0, 1, 0, 1, 0. 

Rj^ ( t) = 32, 1, 0, 1, 0, -3, 0, 5, 0, 1, 0, -9, 0, 

-1 , 0,-1, 16, +1 , 0 , +1 , 0 , +1 , 0 , +1 , 0 , 3 , 
0, -5, 0, -1, 0, -1, 0 
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and sum 



S = Rq(t) + Rq(t) = 64, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, 



0, 


0, 


0, 


0, 


0, 


32, 0, 0 


0, 


0, 


0, 


0, 


0, 


0 5 0 5 0 5 



0 , 0 



with main lobe level 64 and positive sidelobe levels 32 at 
T = - 16 as predicted by the rule. 

For this example, the IBM- 360 computer was used. 



2 . 



^ 2 

Sidelobes - K Away 



Let A, B be complementary sequences of length n. Let 

be a new sequence generated by taking only the y first digits 

of A (truncating A after its y first digits) . Also let $ be a 

new sequence generated by taking only the j fi^st digits of B 

(truncating B after its y first digits). Also let 4^, f be 

complementary of length Then the new sequences of length 

n 3n 

K = n + y = — are 

C = a4^ 

D = B$ 

and C and D are almost complementary sequences. The sum of 

their autocorrelation functions has two positive sidelobes 

1 . +2 
only of magnitude that of the main lobe and at t = - ^ 

In this case there is an improvement relative to the previous 

case in that the sidelobes are smaller compared to the main 

lobe and farther removed from the main lobe. 



The sequences A, B of length n can be constructed by using 
two complementary sequences of length ^ a.ccording to the rule 
in Section III.B(a) . Then fhe sequences ^ are always 
complementary. For example, consider the two complementary 
sequences of length n = 8 

A = 00011101 
B = 01001000. 

The new sequences of length 12 are 

C = A^ = 000111010001 
D = B$ = 010010000100 

with autocorrelation functions 

R^(t) = 12, 1, 0, -5, 0, -5, 0, 1, 4, 1, 0, -1, 0. 
Rj^(t) = 12, -1, 0, 5, 0, 5, 0, -1, 4, -1, 0, 1, 0. 

and s um 

I = Rq(t) + Rq(x) = 24, 0, 0, 0, 0, 0, 0, 0, 8, 0, 

0 , 0 , 0 

1 

with main lobe level 24 and sidelobe levels J (24) = 8 at a 
distance x = - 8 . 

This example is illustrated in Fig. 11. Another 
example is 

A = 01001000000111010100100011100010 
B = 01001000000111011011011100011101 



48 



i 




-1 
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#1 



J 

1:1 

J 



I 



I'; 

ti 







A, B are complementary of length n = 32. 

The new sequences of length K = 48 are 

C = = 01001000000111010100100011100010 

0100100000011101 

D = = 10111000111011011011100000010010 

1011100011101101 

with autocorrelation functions 

R(.(t) = 48, 1, 2, 3, 0, -1, -2, 1, 0, 3, -2, 5, 0, 

+9, 2, -5, 0, -5, 2, 9, 0, 5, -2, 3, 0, 1, 

-2, -1, 0, 3, 2, 1, 16, 1, 2, 3, 0, -1, -2, 

1, 0, -1, -2, 1, 0, -3, 2, -1, 0. 

Rj^(t) = 48, -1, -2, -3, 0, 1,. 2, -1, 0, -3, 2, -5, 

0, -9, -2, 5, 0, 5, -2, -9, 0, -5, 2, -3, 

0, -1, 2, 1, 0, -3, -2, -1, 16, -1, -2, -3, 

0, 1, 2, -1, 0, 1, 2, -1, 0, 3, -2, 1, 0 . 

and sum 

I = R(.(t) + R£)(t) = 96, 0, 0, 0, 0, 0, p, 0, 0, 0, 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

32, 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 0 , 0 . 

The main lobe has level 96 and the sidelobes have levels J (96) = 32, 

at a distance ~ “ " 3 " ~ “ predicted by the rule. 
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Here also the IBM- 360 computer was used because of the 
length of the sequence. 

B. NEGATIVE SIDELOBE SEQUENCES 

Here two cases are developed similar to the ones considered 
before . 

+ 1 

1. Sidelobes - "2 K Away 



If A, B are complementary sequences of length n and the 

new sequences C = AA and D = BB of length K = 2n are formed, 

where A, B represent the complements of A, B respectively, then 

the sum of their autocorrelation functions gives two negative 

sidelobes only with magnitudes half that of the main lobe and 

+ 1 ■ 

at a distance x = - K, 

For example 

A = 1011 
B = 1110 

are complementary. Then 

C = AA = 10110100 and 
D = BB = 11100001 

are almost complementary with autocorrelation functions 

R^(t) = 8, -3, 0, 3, -4, 1, 0, -1, 0 
Rj^(x) = 8, 3, 0, -3, -4, -1, 0, 1, 0 
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m 



and sum 



Z = "®’ °> °* 

Here the sidelobes are negative with levels half that of the 
main lobe att =ljK=^j9 =^4. This is illustrated 
in Fig. 12. 

Another example, for which the IBM- 360 computer was 

used is: 



A = 0001001011100010 
B = 0100011110110111 

A, B are complementary of length n = 16 . Constructing 
C = AA, D = BB gives 

C = 00010010111000101110110100011101 
D = 01000111101101111011100001001000 

Sequences C and D are almost complementary of length 
K = 32 with autocorrelation functions 

R^(t) = 32, -3, 0, -3, 0, -7, 0, 1, 0, -5, 0, 11, 0, 
3, 0, 3, -16, 1, 0, 1, 0, 1, 0, 1, 0, 3, 0, 
-5, 0, -1, 0, -1, 0 

R^(t) = 32, 3, 0, 3, 0, 7, 0, -1, 0, 5, 0, -11, 

0, -3, 0, -3, -16, -1, 0, -1, 0, -1, 0, 

-1, 0, -3, 0, 5, 0, 1, 0, 1, 0 
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Fig. 12. Negative sidelobes at x 
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and Sinn 



Z = Rq(t) + Rj)(t) - 64, 0, 0, 0, 0, 0, 0, 0, 0, 0, 



0, 


0, 


0, 


0, 


0, 


0, -32, 


0, 


0, 


0, 


0, 


0, 


0, 0, 0 



0 , 0 



with main lobe level 64 and negative sidelobe levels 32 at 
T = - 16 as predicted by the rule . 

+ 2 

2 . Sidelobes - K Away 



Let A, B be complementary sequences of length n. Let 

$ be a new sequence generated by taking only the complement of 

n . . . 

the ^ first digits of A. 

Let f be a new sequence generated by taking only the 

complement of the j first digits of B. Also, let f be 

n. - 

complementary of length 2> Then the new sequences C = A^ and 

* ^ 3n 

D = Bf of length K = n + 'f = — will be almost complementary 

and the sum of their autocorrelation functions will have two 

negative sidelobes only with magnitudes one-third that of the 

+ 2 

main lobe at x = - j R- 
For example, 

A = 11101101 
B = 10111000 

are complementary sequences of length n = 8 . The new sequences 
C = A$ = 111011010001 and D = B| = 101110000100 of length K = 12 
are almost complementary with autocorrelations 
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R(.(t) = 12, -1, 0, 1, 0, -1, 0, 1, -4, -1, 0, 1, 0 
Rj^Ct) = 12, 1, 0, -1, 0, 1, 0, -1, -4, 1, 0, -1, 0 

and sum 



^ ~ ~ 

0 , 0 , 0 

The only sidelobes are negative with level (24) = 8 at 

^12) = - 8. This is illustrated in Fig. 13. 

Another example obtained with the use of the IBM- 360 
computer is 



A = 01001000000111010100100011100010 
B = 01001000000111011011011100011101 

A, B are complementary of length n = 32 . The new se- 
quences of length K = 48 are 

A = A$ = 0100100000011101010010001110001010110 
11111100010 

B = B| = 1011100011101101101110000001001001000 
11100010010 . 

The sequences C and D are almost complementary with auto- 
correlations 



Rq (t) — 48, — 1 , 6, — 3, 0, 1, — 6, — 1 , 0, 1, -6, 1, 

0 , 3 , 6 , 1 , 0 , — 1 , — 6 , — 3 , 0 , 1 , 6 , — 1 , 0 , 

1, 6, -1, 0, 3, -6, 1, -16, -1, -2, -3, 0, 

1, 2, -1, 0, 1, 2, -1, 0, 3, -2, 1, 0 
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1 

-4 





Fig. 13. Negative sidelobes at x = - j K. 
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Rp(x) = 48, 



1, -6, 3, 0, -1, 6, 1, 0, -1, 6, 1, 0, 



-3, -6, -1, 0, 1, 6, 3, 0, -1, -6, 1, 0, 
-1, -6, 1, 0, -3, 6, -1, -16, 1, 2, 3, 0, 
-1, -2, 1, 0, -1, -2, 1, 0, -3, 2, -1, 0 

and s\im 



2 = Rq(t) + Rp(T) = 96 


o 

r^ 


, 0, 0, 0, 0 


, 0, 0, 0, 0 


0, 


0, 


0, 0, Oj 0, 


0, 0, 0, 0, 


0, 


0, 


o 

o 

o 

o 


0, 0, 0, 0, 


0, 


0, 


-32, 0, 0, 


0, 0, 0, 0, 


0, 


0, 


0, 0, 0, 0, 


0 , 0 , 0 , 0 



with a main lobe level of 96 and negative sidelobe levels of 
32 at T = - 32 as predicted by the rule. 

In signal detection applications almost complementary 
sequences with negative sidelobes offer better noise immunity 
than those with positive sidelobes. The negative sidelobes 
can be removed with an envelope detector. 
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V. NUMBER OF DIFFERENT AUTOCORRELATION FUNCTIONS IN ALL 



BINARY SEQUENCES OF FIXED LENGTH 



Binary sequences useful for communications or ranging 
purposes have autocorrelation functions with small sidelobes . 
Small can be defined in terms of a predetermined level. 

How can these sequences be found? At present, the only 
way to find these sequences is to form the autocorrelation 
functions of all possible sequences of a given length and 
then select the desirable ones. This is a tedious task 
specially for long sequences, since there are 2^ possible 
binary sequences of length n. However, as shown in this 
section, many of these 2^ sequences have the same auto- 
correlation function. 

In general there are four sequences of given length which 
have the same autocorrelation function: 

1. The sequence itself. 

2. The sequence obtained by reversing the original 
sequence . 

3. The sequence obtained by complementing the original 
sequence . 

4. The sequence obtained by complementing and reversing 
the original sequence. 

For example, consider the sequence of length n - 4: 

A = 1101 
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Reversing A gives 



= 1011 = B. 

If A is complemented, results in 

A = 0010 = C. 

If A is reversed and complemented, gives 

A = 0100 = D. 
r 

All these sequences A, B, C, and D have the same auto- 
correlation function 

These results are easy to prove. Consider each case 
separately . 

1. When a sequence is reversed and its autocorrela- 
tion function taken, this is exactly the same as if the auto- 
correlation function of the original sequence was taken, 
since the formation of the autocorrelation function can be 
considered as being accomplished by "sliding" the sequence 
past itself either to the right or to the left. So, "sliding" 
to the right for A is equivalent to "sliding" to the left 

for A^. 

2. The autocorrelation function of a sequence is 
generated by forming products and adding them. The general 
form of one of these products is 
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a . a . 

1 1 



a^ and a^ can have the values 



a^ = 1 or -1 



a^ = 1 or -1, 



Here the values 1 and -1 are used since multiplication is 
considered in the formation of the autocorrelation function 
(if modulo 2 addition were considered, the values 1 and 0 
would be used) . 

Now the possible values of the product a^ a^ are 



a. 

1 


a . 

3 


= 


(i)(i) = 1 




or 


a . 

1 


a . 

3 


= (i)(-i) 


= -1 


or 


a . 

1 


a - 

3 


= (-D(l) 


= -1 


or 


a . 

1 




= (-D(-l) 


= 1 



If the sequence is complemented, 1 is replaced by -1, and 
-1 by 1. So the possible values of the product a^ a^ are 
respectively 



a . 


a . 

3 


= (-D(-l) = 1 


1 




or 


a . 


a^ = (1)(-1) = -1 




1 


or 


a . 


a. = (1)(-1) = -1 




1 


] 


or 


a . 


a. = (1)(1) = 1. 




1 


D 



The values of the product 3 ^ a.^ are the same as before 
for all i and j and so, the autocorrelation function remains 
the same . 
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3 . Since the autocorrelation function is the same if 
the sequence is reversed or complemented, it follows that it 
will be also the same if the sequence is reversed and comple- 
mented . 

Some sequences are their own reverse. These sequences 
are called here symmetric and will be denoted by the letter S . 

For example the sequence A = 1001, when reversed gives 
= 1001 = A. In this case only the sequence itself and its 
complement need to be considered. The sequences A = 1001 and 
A = 0110 have the same autocorrelation function. 

For some sequences, complementing and reversing gives the 
same sequence. These sequences are called here R sequences. 

For example, the sequence A = 000111 when reversed and 
complemented gives A^ = 000111 = A. In this case only the 
sequence itself and its reverse need to be considered. The 
sequences 000111 and 111000 have the same autocorrelation 
function . 

It is assumed that the number of sequences having the 
same autocorrelation function is either two or four. It has 
not been shown though that there does not exist any other 
number of sequences such as 3 or 5 or 6, etc., that have the 
same autocorrelation function. In this work, it has been 
verified by exhaustive computer search for sequences of various 
lengths, that the number of sequences having the same auto- 
correlation function is either two or four. 

When an R sequence is reversed and complemented, the ones 
of the R sequence become the zeros of the reverse complement . 
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So, an R sequence has always the same number of ones and zeros. 
Therefore, an R sequence is always of even length. 

For example, for length n = 3, there will not be any R 
sequence, since n is odd. 

So far it has been established that: 

(1) There are only four sequences with the same auto- 
correlation function, provided these sequences are not S or R. 

(2) There are only two sequences with the same auto- 
correlation function, if these sequences are R or S . 

(3) S sequences can be of any length. 

(4) R sequences can be only of even length. 

For example, all the sequences of length n = 2 and their 
autocorrelation functions are: 



Sequence Autocorrelations 



00 


2 , 


1 , 0 


01 


2 , 


- 1 , 0 


10 


2 , 


- 1 , 0 


11 


2 , 


1 , 0 



Here there are two S sequences, 00 and 11, with the same 
autocorrelation function = 2, 1, 0. Also, there are two R 
sequences, 01 and 10, with the same autocorrelation function 
= 2 , - 1 , 0 . 

Another example of length n = 3 is : 
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Sequence 



Autocorrelation 



000 


3, 


2, 1, 0 


001 


3, 


0, -1, 0 


010 


3, 


-2, 1, 0 


oil 


3, 


0, -1, 0 


100 


3, 


0, -1, 0 


101 


3, 


-2, 1, 0 


110 


3, 


0, -1, 0 


111 


3, 


2, 1, 0 



Here there are four sequences which are neither R nor S, 
001, Oil, 100. 110 with the same autocorrelation function 
3, 0, -1, 0. Also there are four S sequences: 000, 111, 010, 

101. The sequences 000 and 111 have the same autocorrelation 
function 3, 2, 1, 0 and the sequences 010 and 101 have the 
same autocorrelation function 3, -2, 1, 0. In this example 
there is no R sequence since n = 3 is odd. 

Next the exact number of R and S sequences will be 
established among all the possible sequences of length n. 

The maximum number of binary sequences of length n is given 
by 2^. The sequences of n = 1 are 0 and 1. The number of 
S sequences is two. So S = 2 here;(0, 1). For n = 2 all 
sequences are 00, 01, 10, 11. Here S = 2: (00, 11), R = 2: 

(01, 10), and S + R = 4. For n = 3, the possible sequences are 
000, 001, 010, oil, 100, 101, 110, 111. Here S = 4:(000, 010, 
101, 111), R = 0 and S + R = 4. Going to n = 4 , it can be 
seen that S = 4 , R = 4 and S + R = 8 . 
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The relation between n and S and R can be derived by 
considering the mechanism of moving from a sequence of length 
n to the next one of length n + 1. 

For example, the sequences of n = 2 are formed by taking 
the sequences of length n = 1 and adding in front of each of 
them a zero and a one, one at a time, so the number of sequences 
for n = 1 is doubled and all the possible sequences of length 
n = 2 are formed. By doing so, the following can be noted: 

(a) When moving from n even to n+1 which is odd, the 

S + R sequences in n is equal to the S sequences in n+1. 

(b) When moving from n odd to n+1 which is even, the 

S in n is half the S + R in n+1. 

For example, for n = 4 there are S + R = 8 sequences and 
for n = 5 there are S = 8 , but moving to n = 6 gives S + R = 16 . 

So, moving from an even length to the next keeps the number 
S + R, but moving from an odd length to the next doubles the 
S + R number. 

Using formulas it can be written: 

(a) n = even, (S + R)^ = 

(b) n = odd, = y (S + 

Now, considering the following short table giving the S + R 
terms of all possible sequences up to length n - 10 , the 
relation between n and S + R can be obtained. 
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I 



$ -m 




n 



S + R 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

For n even, S + R 
n odd, S + R = 



2 

4 

4 

8 

8 

16 

16 

32 

32 

64 

n + 2 

2 ^ and for 
n+1 



The number of different autocorrelation functions contained 

in all possible sequences of length n is derived as follows. 

If this number is denoted by T, two cases are considered. 

(a) n is even. Here all the possible sequences number 

n+2 

n 2 

2 . Also the number of S and R sequences is S + R = 2 . By 

deducting S + R from 2^, a number of sequences equal to 
n + 2 



Tl 9 

2^-2 is obtained. It was established before that since 

n+2 



these remaining 2^ - 2 ^ sequences contain no R or S sequences, 

n+2 

2 ^ - 2 ^ • • 
they have jr different autocorrelation functions. The 

2 2 ^ 

2 S and R sequences have — 2 — different autocorrelation 
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functions. So finally, all the 2^ sequences have 



n+2 



n+2 



2n _ 2 2 



“ + — 2 — different correlation functions. By re- 
arranging this result, the final formula for n even is derived 

(T) = 2 ^~'^ + 2 ^^^ 

even 



(b) 



n is odd. Similarly all possible sequences here 



n+1 



n+1 



number 2^. Also, S + R = 2 ^ , Deducting the 2 ^ S + R 



n+1 

n n 2 

sequences from 2^ gives 2^^ - 2 sequences with 



n+1 

2n _ 2 2 



different autocorrelation functions. So, all the 2^ sequences 
n+1 n+1 

«n _ 2 2 2 2 

have ^2-r; + — s — different autocorrelation functions 



4 2 

or in the final form 



n-1 



(T) = 2^”^ + - 

^^''odd ^ 2 



Two examples are taken 
(a) n = 7 



<T>odd = 2® " V = 36. 



So the 2^ = 128 possible sequences of length 7 give only 36 
different autocorrelation functions . 



(b) n = 10 
(T) 



= 2^ + ^ = 272 . 
even 2 
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So, by using sequences 10 bits long, at most 272 different 
autocorrelation functions can be obtained from the 2^^ = 1024 
possible sequences. 
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VI . COMPUTER SEARCH FOR "GOOD CODES 



In "this section the computer programs for obtaining sequences 
with "small" autocorrelation sidelobes are discussed. These 
sequences or "good" codes are found by applying mainly the 
results of the previous section. Because of the large number 
of different sequences of length n for even modest values of n , 
it is necessary to use a digital computer to search for "good" 
codes. Practically, an exhaustive search is limited to n ~ 20 
with present digital computers. 

For small lengths (up to n = 10) it is possible that a 
programmable calculator can be used to find the autocorrelation 
function of one sequence at a time. This was done with a TI-59 
programmable hand calculator by storing each element of the 
sequence in a memory location. Then the autocorrelation function 
was formed by multiplication of the proper elements in each 
position. In this case, +1 and -1 is used for the elements 
of the sequence. 

The algorithm for computing the autocorrelation function 
of a sequence on a large computer (IBM- 360) is constructed. 

This is given in Program 1 on page 111. In this program the 
autocorrelation function of only one sequence can be computed. 

The sequence has to be punched on a computer card. A small 
modification gives Program 2 on page 112 which gives 
the option to find the autocorrelation functions of any number 
of sequences. Each of these sequences has to be punched on 
a separate card. 
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In these computer programs and also all the next ones, the 
sequences are represented with zeros and ones. The algorithm 
for finding the autocorrelation function has been made by 
comparing the number of like and unlike elements in every 
position . 

To avoid punching the sequences on the cards, a program is 
created to generate automatically all the possible 2^ sequences 
at length n. This is accomplished by counting in binary from 
0 to 2^ and thus generating all the binary numbers of length n. 
This program is combined with the program for the computation 
of the autocorrelation function. So, every time a sequence is 
generated, its autocorrelation function is formed. This is 
Program 3 on page 113. 

Since only the different autocorrelation functions are of 
interest here, a program is written to select only those codes 
having different autocorrelation functions. To understand the 
operation of this program, consider the following example. 

Take all the sequences of length n = 4 and the corresponding 
autocorrelation functions . 
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Sequence Autocorrelation 



0000 


4, 


3, 2, 1, 0 


0001 


4, 


1, 0, -1, 0 


0010 


4, 


-1, 0, 1, 0 


0011 


4, 


1, -2, -1, 0 


0100 


4, 


-1, 0, 1, 0 


0101 


4, 


-3, 2, -1, 0 


0110 


4, 


-1, -2, 1, 0 


0111 


4, 


1, 0, -1, 0 


1000 


4, 


o 

1 — 1 

1 

o 


1001 


4, 


o 
1 — 1 

CM 

1 

rH 

1 


1010 


4, 


-3, 2, -1, 0 


1011 


4, 


O 
« — 1 
o 

1 — 1 

1 


1100 


4, 


1, -2, -1, 0 


1101 


4, 


-1, 0, 1, 0 


1110 


4, 


1, 0, -1, 0 


1111 


4, 


3, 2, 1, 0 



4 

There are a total of 2 =16 sequences . After the eighth 

sequence 0111, the other sequences are complements of the 
first eight ones. So, they give no new autocorrelation 
function, and therefore they do not need to be taken under 
consideration. In the first eight sequences there are two 
pairs with the same autocorrelation function 0001 and 0111 
and also 0010 and 0100. So, only six different autocorrela- 
tion functions remain. The general way to proceed is to take 

2 ^ 

each sequence in the first — , reverse it, complement it and 
reverse complement it and then keep only the original sequence 
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and reject the others, since they have the same autocorrela- 
tion function. Here two cases are considered. 

(a) The sequence ends in 0. 

In that case its complement and reverse comple- 
ment will start with one, so they belong to the sequences after 
2 ^ 

the first — ones and need not to be considered. In that case 
only the reverse of the sequence is taken. 

(b) The sequence ends in 1. 

In that case its reverse and complement will start 

2 ^ 

with one, so they belong to the sequences after the first — ones 
and need not to be considered. In that case only the reverse 
complement is taken. 

Now, the program takes each sequence as it is generated and 

reverses it if it ends in 0 , or reverse complements if it ends 

in 1. If the resulting sequence represents a smaller binary 

number than the original one, this means that the resulting 

sequence was generated before and its autocorrelation function 

already taken, so there is no need to be taken again, and the 

program goes to the next sequence. The procedure is repeated 

2 ^ 

until all the first — sequences are finished. This way only 
the codes having different autocorrelation functions are listed 
in the printed output. 

For example, in the case of n = 4 the result is as follows: 
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Sequence Autocorrelation 



0000 


4, 


3, 2, 1, 0 


0001 


4, 


1, 0, -1, 0 


0010 


4, 


-1, 0, 1, 0 


0011 


4, 


1, -2, -1, 0 


0101 


4, 


-3, 2, -1, 0 


0110 


4, 


-1, -2, 1, 0 



And the number of different autocorrelation functions is 6 , 
which is in agreement with the formula 



(A) 



even 



2 



n- 



2 



+ 



2n/2 




6 



The program which generates automatically the sequences 
of length n and computes only the different autocorrelation 
functions is Program 4 on page 114 . An example for n = 10 is 
given on page 84. 

All the different autocorrelation functions are not needed. 
Only those with small sidelobes . So a filtering procedure has 
to be introduced in the program to keep only those auto- 
correlations which have sidelobes equal or smaller than a pre- 
determined level. Program 5 on page 115 generates automatically 
the sequences of length n and prints only those different auto- 
correlation functions with sidelobe levels equal or less than 1. 
For this case the lengths n = 10 , 11, 12, 13, 15, 17 and 20 
were examined, and the number of different autocorrelation 
functions found with sidelobe levels equal or less than 1 were, 
respectively, 11, 1, 16, 31, 0, 40, 3. 
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It is interesting to note that for n=ll there is only one 
such autocorrelation which is a Barker code (a sequence with 
sidelobes between +1 and -1). For n=15 no such codes exist, 
and for n=20 only three were found. 

For the case n=20 it was not possible to make an exhaustive 
search because of the computer time required. Using 30 
minutes of computer time, only three such autocorrelation 
functions were found. It is estimated that about 2 hours of 
computer time is required to make an exhaustive search. 

It is evident that for lengths greater than n=20 even with 
a large computer an exhaustive search is impractical. 

There is though a way to search regions of big sequences . 
Computer Program 6 on page 116 was used to search for a region 
of the n=20 case with no new results. 

The results for the cases n=10, 11, 12, 13, 17, 20 
are listed on pages 96, 98, 99, 101, 105, 110, respectively. 
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VII. APPLICATIONS 



Complementary sequences, almost complementary sequences 
and codes with small sidelobes can be used in communications, 
ranging and spread spectrum systems . This section lists some 
possible applications. 

Two complementary sequences can be transmitted simultane- 
ously using a quadriphase phase-shift keying (QPSK), for 
example . Two complementary sequences A and B are applied to 
a QPSK modulator. The output V^(t) is a sine wave with a 
phase which can have one of four values . In the receiving 
system, V^(t) is demodulated. The outputs of the demodulator 
are the sequences A and B. Each sequence is processed by a 
matched filter and the outputs of the two matched filters 
added. Then the output of the summer will have one main lobe 
and no sidelobes. Fig- 1*+ illustrates this system. 

Almost complementary sequences with positive sidelobes 

can be used to measure doppler. In Fig. 15 such a scheme is 

shown. The almost complementary sequences A and B of length K 

+ 2 

give after summing positive sidelobes at t = - -^ K from the 
main lobe. Each of them is processed by a matched filter and 
the two outputs of the matched filters added. The output of 
the summer is connected to a counter in such a way that the 
first pulse (main lobe) enables the counter and the second 
pulse (sidelobe) inhibits the counter. When doppler is 



75 



A 



B 



' QPSK 


V^(t) 


Modulator 





(a) Transmitter 



V^(t) 


QPSK 


^ 


M. F. 

A 


^ ^ 


Demodulator 










DO 


M. F. 

EL_ 




no sidelobes 



(b) Receiver 
Fig. 14. QPSK System. 
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Ra(t) 



+ Rj(T) 



(b) Receiving system. 

Fig. 15. System to measure doppler. 
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introduced to the system, the distance between the main lobe 
and the sidelobe will change by an amount proportional to the 
doppler. The contents of the counter is, then, a measure of 
the doppler. 

Codes with small autocorrelation sidelobes can be used in 
spread spectrum systems. They can be used as "chip" sequences 
and as means for synchronizing the remote oscillators in such 
systems. They can also be used in surveillance and ranging 
systems because they provide accurate and unambiguous time 
measurements . 
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VIII. RESULTS AND CONCLUSION 



The study presented in this thesis had as a main objective 
the discovery of means of generating "good" codes . Good codes 
imply binary sequences having autocorrelation functions with 
small sidelobes or no sidelobes at all. 

The results of the search for good codes in this study are 
the following: 

(1) Almost complementary sequences were discovered. 

These are constructed by using complementary sequences . The 
sum of the autocorrelation functions has two sidelobes only of 
predictable level, polarity and position. 

(2) Computer programs for obtaining sequences with 
autocorrelation functions having sidelobe levels less than 
predetermined desirable levels were prepared and used. These 
computer programs reveal that for sequences of length n = 15, 
there are no codes with sidelobes less than or equal to 1. 

There is only one such code for n = 11 (which is a Barker code) 
and for n = 20 three such codes were found without exhausting 
all possibilities. For lengths greater than n = 20, an 
exhaustive computer search is impractical. A partial search 
can be made though by searching regions of these sequences 
with a special computer program. 

The search for good codes resulted also in the following: 

(1) A formula which gives the number of the different 
autocorrelation functions for all possible sequences of length n. 
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(2) A computer program which can generate automatically 
all sequences of length n and give their different autocorrela- 
tion functions. 

Some application of these codes were also considered. 

There are some suggestions for further research. 

(1) Supplementary and cyclic complementary sequences 
and their properties could be further investigated and appli- 
cations developed. 

(2) Applications of the sequences presented here 
could be implemented with hardware. 
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APPENDIX A 



COMPUTER OUTPUTS AND PROGRAMS 



Six computer programs are included in this report and 
seven computer outputs which are the results of the search 
for the "good" codes. In the programs, the sequence length 
is denoted by L. 

Program 1 computes the autocorrelation function of one 
sequence, by comparison of the elements of the original 
sequence and its shifted replica. The length L of the sequence 
and the sequence (CODE (I)) have to be punched on separate 
cards . 

Program 2 computes the autocorrelation function of any 
number of sequences of the same length. The number of the 
sequences has to be specified in the loop DO 150. The length 
of the sequences and each of the different sequences have to 
be punched on separate cards . 

Program 3 generates all the possible sequences of length 
L and computes their autocorrelation function. Only the 
length L has to be punched on a card. Everything else is 
done automatically. Thus, punching of the sequences on cards 
is avoided. 

Program 4 computes only the different autocorrelation 
functions of all sequences of length L. This is done by the 
algorithm explained in Section VI. The program produces 
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everything automatically and only the length L has to be 
specified and punched. 

Program 5 computes the different autocorrelation functions 
of all sequences of length L which have sidelobes less or 
equal to one. The length L has to be specified. 

Program 6 searches only a region of all sequences of 
length L and computes the different autocorrelation functions 
with sidelobes less or equal to one which exist in this region. 
Here the length L and the starting sequence (SCODE (I)) have 
to be specified and punched on different cards each. All the 
sequences before the starting sequence are ignored. This was 
done because for long sequences, a complete computer search 
is impractical. 

There is a computer output which lists the different auto- 
correlation functions for all sequences of length 10 . The 
rest of the outputs list the "good" codes and their autocorrela- 
tion functions for lengths 10, 11, 12, 13, 17 and 20. The codes 
listed in all computer outputs are printed with the most signif- 
icant bit on the right . 
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PRINCIPAL VARIABLES USED 



CODE 

COR 

L 

ICODE 

JCODE 

SCODE 



binary sequence 

autocorrelation function of a binary sequence 
length of a binary sequence 

variable used to complement a binary sequence 
variable used to reverse a binary sequence 
starting sequence when a region of the binary 
sequence of length L is searched 
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PROGRAM «1 



COMPUTATIOM OF THE AUTOCORRELATION FUNCTION 
OF CNE BINARY SEQUENCE 

INTEGER COOE( 50) ,C0R(50 ) 

READ(5,10)L 
10 FORMAT! 12) 

REAO(5,1000 ) (CODE! I), 1=1, 48) 

1000 FORMAT! 5011) 

DO 40 1=1, L 
N=0 

DO 30 J = I ,L 

IF!CODE! J ).E0.C00E!J-! I-l) ) ) GO TO 20 
N=N-1 
GO TO 3 0 
20 N=N+1 
30 CONTINUE 
40 COR! I) =N 
COR !L+1 ) = 0 
K=L+1 

WRITE!6,100) (CODE! I) ,I =l ,L) 

WRITE!6,100) !COR(I ), I=1,K) 

100 FORMAT!// ,’ ', 30(13, IX)) 

END 
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PROGRAM 4Z 

COMPUTATION OF THE AUTOCORRELATION FUNCTION 
OF y BINARY SEQUENCES 

INTEGER COOE( 50) ,C0R(50) 

REA0(5 ,10 )L 
10 F0RMAT(I2) 

00 150 M= I, A 

REAO(5tlOOO) (C006( I), I=1,A8» 

1000 F0RMAT(50I1) 

DO AO 1=1, L 
N=0 

DO 30 J=I ,L 

IF(CODE(J i.EQ.CODEI J-{ I-in ) GO TO 20 
N=N-1 
GO TO 30 
20 N=N + 1 
30 CONTINUE 
AO CORU )=N 
COR(L+l ) = 0 
K=L+1 

WRITE(6,100) (CODE( I) ,I =l ,L) 

WRITE(6,100) (COR(I ) ,I = 1,K) 

100 FORMAT!//,' ’,30(13, IX)) 

150 CONTINUE 
END 
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PROGRAM *3 



AUTOMATIC PRODUCTION OF ALL SEQUENCES 
OF LENGTH L AND COMPUTATION OF 
THEIR AUTOCORRELATION FUNCTIONS 

INTEGER C0DE(20)t COR (20 ) 

READ(5»10)L 
10 F0RNAT(I2) 

NN=2**L 

GENERATE SEQUENCES 

00 100 H=ltNN 
DO 20 1=1,20 
20 C00E(H=0 
NI=II-1 
DO 30 M=1,L 
JJ=L-M 

IF (NI.LT. 2**JJ) GO TO 30 
CODE( JJ+1 ) = 1 
NI=NI-2**JJ 
30 CONTINUE 

COMPUTE AUTOCORRELATIONS 

00 60 I=1,L 
N=0 

DO 50 J=I ,L 

IF(CQDE(J ).£Q .CODE( J-(I-l)) ) GO TO AO 
N=N-1 
GO TO 50 
A-0 N=N+1 
50 CONTINUE 
60 COR(I )=N 
C0R(L+1 )=0 
K=L+1 

WRITE(6,200) (CODE ( I), 1=1, L) 
WRITE(6,200) (COR( I ) , 1=1 ,K) 

100 CONTINUE 

200 FORMAT!//,' ',20(I2,2X)) 

END 
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I 








PROGRAM 



10 



AUTOMATIC CALCULATION OF THE DIFFERENT 
AUTOCORRELATION FUNCTIONS IN ALL POSSIBLE 
SEQUENCES OF LENGTH L 

INTEGER C00E(20),C0R(21)f IC00E(20) ♦JC00E(20) 
READ(5 ,10 )L 
FORMAT! 12 ) 

NN=2**<L-1) 

GENERATE SEQUENCES 



DO 100 11=1, NN 
DO 20 1=1,20 
20 C0DE(I)=0 
M =11-1 
DO 30 M=l,L 
JJ=L-M 

IF(M.LT.2*»JJ) 
CODEUJ+l )=1 
NI=NI-2**JJ 
30 CONTINUE 



GO TO 30 



TEST IF SEQUENCE ENDS IN 0 OR 1 

DO 31 M=1.L 

31 ICODEIM )=CODE (M ) 

IF (COOEd I .EQ.Ol GO TO 33 

COMPLEMENT SEQUENCE 

DO 32 M=1 ,L 

IC0DE(M)=IC0DE{M)+1 

IF (ICODEI M) .GT.l ) IC0DE(M) = 0 

32 CONTINUE 

REVERSE SEQUENCE 

33 DO 34 M=1 ,L . 

34 JCCDE(L-M+1 )=ICODE(MJ 

IF( JCOOEI L-M+1)-C0DE(L-M+1) ) 100, 35,36 

35 CONTINUE 

36 CONTINUE 

COMPUTE AUTOCORRELATION 

DO 60 1=1 ,L 
N=0 

IF (CODEU )*^EQ .CODE ( J-n-1) ) ) GO TO 40 
N=N-1 
GO TO 50 
40 N=N4l 
50 CONTINUE 
60 COR(I)=N 
C0R(L+1 )=0 

70 WRITE(6,200) (CODE ( I ) , I =1 ,L ) , (COR ( I ) ,I =1 , K) 
100 CONTI NUE . . 

200 FORMAT!/,’ ',41(12, IX)) 

END 
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I 1 



PROGRAM 



10 



20 



30 



31 



32 



COMPUTATION OF THE DIFFERENT 
AUTOCORRELATION FUNCTIONS 
OF ALL SEQUENCES OF LENGTH L 
WITH SIDE LOBES LESS OR EQUAL TO 1 

INT EGER CODE! 20 ), COR ( 21 ) , ICOOEI 20 ) , 

REAO(5,10)L 

FORMAT! 12 ) 

NN=2**(L-1 ) 

GENERATE SEQUENCES 

DO 100 11=1 ,NN 
00 20 I=lf20 
CODE! I)=0 
NI=II-1 
00 30 M=1»L 
JJ=L-M 

IF (NI.LT.2’«'»JJi GO TO 30 
CODE! JJ>1 ) = 1 
N1 = N I-2**JJ 
CONTINUE 

TEST IF SEQUENCE ENOS IN 0 OR 1 



00 31 M=1,L 
ICODE!M ) = COOE!MI 
IF (COOEIl ) .EQ.OI 



GO TO 33 



COMPLEMENT SEQUENCE 



00 32 M=1,L 
ICOOEIM )=IC0DE!MJ+1 
IF ! ICODECM) .GT.l ) IC00E(MI = 0 
CONTINUE 

REVERSE SEQUENCE 



33 

34 



35 

36 



00 34 M=1»L 
JCQDE(L-M+1 )=ICOOE!M) 

00 35 M=1 ,L 

IF! JCODE !L-M + 1 )-CODE (L-M+1 ) ) 100, 35. 

CONTINUE 

CONTINUE 

COMPUTE AUTOCORRELATION 



00 60 1=1, L 
N=0 

DO 50 J=I ,L 

IF(CODE!J I.EO.CODE! J-U-D) ) GO TO 
N=N-1 
GO TO 50 
40 N=N+1 
50 CONTINUE 

KEEP SEQUENCES WITH 
SI0EL08ES LESS OR EQUAL TO 1 



55 

60 



70 

100 

200 



IF !N.EQ.L ) GO 
IF! N.GT.l ) GO 
COR !I)=N 
COR ! I)=N 
COR !L+1 J =0 
K = L+1 

WRITE!6,200) 
CONTINUE 
FORMAT!/,' ' 
END 



TO 

TO 



55 

100 



!CODE! I),I=1,U, !COR( 
Al! I2,1X) ) 



JCOOE! 20) 



36 



AO 



I ), I=1,K ) 
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PROGRAM «6 



COMPUTATION OF DIFFERENT 
AJTOCORRE LATI CN FUNCTI ONS FOR A 
REGION OF THE SEQUENCES OF LENGTH L 
WITH SIDELOBES LESS OR EQUAL TO I 

THE SEARCH STARTS AT THE 
SEQUENCE SCODE(I) , I=1,L 

INTEGER CODE (20) , COR (21) ,ICODE( 20 ) , JCOOE (20) , 
REAO( 5,10 )L 
10 F0RMAT(I2) 

READ STARTING SEQUENCE 



11 

12 



20 



30 



31 



READ(5rll) (SCODE ( I ) ,I =1, L) 
FORMAT (2011 I 
ISTART=0 
DO 12 1=1 »L 

ISTART=ISTART+SCODE(I )*2*»( I-l) 
NN=2**( L-D-l 

GENERATE SEQUENCES 

DO 100 II=ISTART,NN 
DO 20 1=1 ,20 
COOE( I ) =0 
NI = II 

DO 30 M=l,L 
JJ=L-M 

IF (NI .LT .2**J J) GO TO 30 
CODE( JJ+1 )=1 
NI=NI-2**JJ 
CONTINUE 

TEST IF SEQUENCE ENDS I N 0 OR 1 

DO 31 M=1 ,L 
ICOOE(M) = (:ODE (M) 

IF( CODE( 1 ) .EQ .0) GO TO 33 

COMPLEMENT SEQUENCE 



DO 32 M=1 ,L 

IC0DE(M)=IC0DE(M)+1 

IF( ICODE(M) .GT«1) ICODE(M) = 0 

32 CONTINUE 

REVERSE SEQUENCE 

33 DO 34 M=1 ,L 

34 JC0DE(L-M-J-1)=IC0DE( M) 

DO 35 M=1,L 

IF ( JCOOE ( L-M+D-COOE (L-M+1 ) )1 00,35 ,35 

35 CONTINUE 

36 CONTINUE 

COMPUTE AUTOCORRELATION 



DO 60 
N=0 



1=1 ,L 



IMCODEU fi'E0.C0DE( J-(I-D) ) GO TO 40 



40 

50 



N=N-1 
GO TO 50 
N=N + 1 
CONTINUE 

KEEP AUTOCORRELATIONS 
WITF SIDELOBES LESS OR 



EQUAL TO 1 



SC00£(20) 
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55 

60 

70 

100 

200 



IF(N .£0 .L ) GO TO 55 

IF ( N.GT.l ) GO TO 100 

COR ( I ) = N 

COR ( I ) = N 

CDR( Ln)=0 

K=L + 1 

WRITE(6t 2 00) (CODE ( I )t I =ltU , (CQR( I),I = 1,K) 
CONTINUE 

FORMAT!/, • ',41(12, IX)) 

END 



117 



LIST OF REFERENCES 



1. Golay, M. J. E., "Multi-Slit Spectometry , " J. Opt. Soc . Am., 
Vol. 39, pp. 437-444, June 1949. 

2. Jauregui, S., Jr., "A Theoretical Study of Complementary 
Binary Code Sequences and Computer Search for New Kernels," 
Ph.D. Thesis, Naval Postgraduate School, Monterey, CA, 

May 1962. 

3. Cooper, G. R., "Probabilistic Methods of Signal and System 
Analysis," pp. 205-210, HPW, 1971. 

4. Ziemer, R. E. and Tranter, W. H., "Principles of Communica- 
tions," pp. 311-314, Houghton Mifflin, 1976. 

5. Golay, M. J. E., "Complementary Series," IRE Transactions 
on Information Theory, Vol. IT-7, pp. 82-87, April 1961. 



118 



INITIAL DISTRIBUTION LIST 



1. Defense Document air ion Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 

3. Department Chairman, Code 62 
Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

4. Associate Professor Glen A. Myers 4 

Code 62Mv 

Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 

5. LT loannis Anastasopoulos H.N. 1 

Panormou 57, Athens 605 

Greece 

6 . Hellenic Navy Headquarters 3 

c/o Embassy of Greece Naval Attache 

2228 Massachusetts Avenue, NW 
Washington, DC 20008 

7. Associate Professor S. Jauregui, Jr. 1 

Code 62Ja 

Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 



No . Copies 
2 



119 



J 



A . f. 1 0 



2 5 5 9 






Thesis 

A4424 

c,l 



179314 

Anastasopoul os 

A theoretical 
study and com- 
puter search 
for binary se- 
quences having 
speci f i c auto-* 
correl ati on 
f uncti ons» 



i Af 



2 5 9 9 



Thesi s 

A4424 

c.l 



Anastasopoul os 
A theoretical 
study and com- 
puter search 
for binary se- 
quences having 
specific auto- 
correl ati on 
functions. 



179314 



thesA4424 

A theoretical study and computer search 




3 2768 000 99121 0 

DUDLEY KNOX LIBRARY 



